Quantum ergodicity on large regular graphs 

^>! Nalini Anantharaman and Etienne Le Masson 

Abstract. We propose a version of the Quantum Ergodicity theorem on large 
QQ ' regular graphs of fixed valency. This is a property of derealization of "most" 

^O _ eigenfunctions. We consider expander graphs with few short cycles (for in- 

stance random large regular graphs). Our method mimics the proof of Quan- 
tum Ergodicity on manifolds : it uses microlocal analysis on regular trees, as 
,-Ch ' introduced in 1231 . 

Oh! 

1. Introduction and main results 

It has been suggested by Kottos and Smilansky that graphs are a good ground 
of exploration of the ideas of "quantum chaos" [201 121) . This means that the spec- 
trum of the laplacian, as well as its eigenfunctions, should exhibit universal features 
that depend only on qualitative geometric properties of the graph. Whereas spec- 
^A ' tral statistics have been extensively studied, both numerically and analytically, the 

■^J" , localization of eigenfunctions have (to our knowledge) only been investigated in a 

,-+ ' few models : the star graphs (both metric and discrete) O [18], the large regu- 

^^ , lar discrete graphs |31l HJ H] , and a family of metric graphs arising from measure 

CO ' preserving 1-dimensional dynamical systems [2]. For the latter, a version of the 

"Quantum Ergodicity theorem" (also known as Shnirelman theorem) has been es- 
tablished. For star graphs, the paper [3] shows on the opposite that "Quantum 
Ergodicity" holds neither in the high frequency limit nor in the large graph limit. 
Furthermore it shows there are eigenfunctions that localise on two bonds of the 
graph. Spectral properties of large regular discrete graphs have been studied in 



> 






, I17|i 133), 130] but eigenfunctions have attracted attention only recently. A 
statistical study of the auto-correlations and the level sets of eigenvectors appeared 
in the papers O 110] that introduce a random wave model (see also |16j for a 
random wave model on metric graphs). The paper [5] has pioneered the study of 
quantum ergodicity on large regular graphs - that is to say, the study of the spatial 
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distribution of eigenfunctions of the laplacian. The result of [5] shows some form 
of derealization of eigenfunctions : 

Theorem 1.1. [5j Let (G n ) be a sequence of (q + l)-regular graphs (with q 
fixed), G n — (V n , E n ) with V n = {1, . . . , n}. Assume tha\^ there exists c > 0, 5 > 
such that, for any k < chin, for any pair of vertices x, y S V n> 

hi 1—s \ 
| {paths of length k in G n from x to y} \ < q \ 2 ' . 

Fix e > 0. Then, if <f> is an eigenfunction of the discrete laplacian on G n and 
if A <zV n is a set such that 

then \A\ > n a — where a > is given as an explicit function of e, S and c. 

A similar form of derealization (but on weaker scales) is established when the 
degree q — q n goes to infinity in [HQ [34] . We also refer to the papers |13L 1141 1121 
111] where various forms of derealization have been established for eigenvectors of 
random Wigner matrices and random band matrices. 

In this paper, our aim is to establish for large regular graphs a result which reads 
like an analogue of the "quantum ergodicity theorem" on manifolds. Compared to 
Theorem II .11 it pertains to a different definition of derealization : derealization is 
now tested by averaging an observable and comparing with the average along the 
uniform measure. As a motivation, let us recall the Quantum Ergodicity theorem 
in its original form. 

Quantum Ergodicity Theorem (Shnirelman theorem, [32|l6l, [36j ). Let 

(M, g) be a compact Riemannian manifoldo Call A the Laplace-Beltrami operator 
on M. Assume that the geodesic flow of M is ergodic with respect to the Liouville 
measure. Let (ipn)nen he an orthonormal basis of L 2 (M,g) made of eigenfunctions 
of the laplacian 

Aip n = -A„Vn, A„ < A„+i — y +oo. 
Let o be a continuous function on M such that J M a(x)d Vol(x) = 0. Then 

(1-1) ^7T7 J^ K^'^nWf)! 2 ^^ 

where the normalizing factor is A^(A) = |{?i,A„ < A}|. Here (ip n ,aip n ) l 2 (m) — 
J M a(x)\^ n {x)\ 2 dVo\(x). 

Remark 1.2. Equation (|l.ip implies that there exists a subset ScH of density 
1 such that 

(1.2) (il>n,aipn) — > / a(a;)dVol(x). 



(Note that (|1.2|l is true even for a function a with non-zero mean.) 

In addition, since the space of continuous functions is separable, one can ac- 
tually find ScNof density 1 such that (TQ|) holds for all a e C°(M). In other 



This assumption holds in particular if the injectivity radius is > clnn. The interest of the 
weaker assumption is that it holds for typical random regular graphs |27| . 
2 We assume that Vol(Af) = 1. 
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words, the sequence of measures (\il>n(x)\ 2 dVol(x)) n es converges weakly to the 
uniform measure dVo\(x). 

Actually, the full statement of the theorem says that there exists a subset ScN 
of density 1 such that 

(1.3) {ll>n,Mn) — ► / <J°(A)dL 

n — >+oo,7ieS J S , M 

for every pseudodifferential operator A of order on M. On the right-hand side, 
a (A) is the principal symbol of A, that is a function on the unit cotangent bun- 
dle S* M , and L is the normalized Liouville measure (uniform measure), arising 
naturally from the symplectic structure of the cotangent bundle. 

Here we consider a sequence of (q+ l)-regular connected graphs (G n )neN, G n = 
(V n , E n ) with vertices V n — {1, . . . ,n} and edge set E n . The valency (g + 1) is fixed. 
We denote by X the (q + l)-regular tree and identify it with its set of vertices, 
equipped with the geodesic distance d%, that we will simply denote by d most of 
the time. We will denote by dc n the geodesic distance on the graph. Each G n is 
seen as a quotient of X by a group of automorphisms T n : G n — T n \X, where we 
assume that the elements of T n act without fixed points. Accordingly, a function 
/ : V n — > C may be seen as a function / : X — > C satisfying f(<y ■ x) = f(x) for 
each x £ X, 7 £ T n . We will denote by D n a subtree of X which is a fundamental 
domain for the action of T n on vertices. 

We consider the stochastic operator acting on r„-invariant functions, 

(1.4) Af{x) = -± ri Y J f(v) 

x~y 

where x ~ y means that x and y are neighbours in the treqj. It is related to the 
discrete laplacian by A — I = A. 

Whereas the Shnirelman theorem deals with the high frequency asymptotics 
(A n — > +°°)j there is no such asymptotic regime for discrete graphs since the 
laplacian is a bounded operator. We will instead work (like in [5]) in the large 
spatial scale regime n — > +00. 

We will assume the following conditions on our sequence of graphs : 

(EXP) The sequence of graphs is a family of expanders. More precisely, there 
exists (3 > such that the spectrum of A on L 2 (G n ) is contained in {1} U [—1 + 
13,1-13] for all n. 

(EUR) For all R, '^ xe "■■ p{ - x > < u — >. where p(x) is the injectivity radius at 
x (meaning the largest p such that the ball B(x,p) in G n is a tree). 

(EIIR) is equivalent to saying that there exists R n — > +00 and a n — > such 

that 

\{x£V n ,p(x) <R n }\ 
< a n . 

n 

In particular, it is satisfied if the injectivity radius goes to infinity (with R n taken 
to be the minimal injectivity radius and a n =0). 

Condition (EXP) replaces the ergodicity assumption in the usual quantum 
ergodicity theorem. 



This is also the (normalized) adjacency matric of the graph G n , but note that this definition 
allows G n to have loops and multiple edges. 
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Example 1. The graph G n can be chosen uniformly at random among the 
(q+ l)-regular graphs with n vertices (see [4j section 2.4 for an introduction to this 
model). We can then take R n = k and na n = AQAkq k for any k — k(n) such that 
kq k n^ 1 — ► 0, and A = A(n) such that A > c > 1 (see [27], Theorem 4). For 

n— >+oo 

instance, we can take k = (1 — (5) log„(n), with < S < 1, and A = 2. In this case 
we have i? n = (1 — S)log q (n) and a„ = 80(1 — S) log q (ri)n~ s . For this choice of 
parameters, (EIIR) is satisfied with a probability tending to 1 when n — > +oo. 
More precisely, this probability is greater than 1 — e~ Cn , for some constant C > 
independent of n. 

Condition (EXP) is also satisfied by these sequences of random graphs : pQ 
proves an equivalence between having a uniform spectral gap and having a uniform 
Cheeger constant. The latter condition was shown to hold generically in [29] . In 
[15] , a spectral gap estimate that is close to optimal is established. 

Example 2. An explicit example of sequence of (<7 + l)-regular graphs to which 
our results apply is given by the construction of Ramanujan graphs of [24 for prime 
q. The sequence obtained satisfies conditions (EXP) and (EIIR) even more strongly 
than the sequences of random graphs of Example [1] A method for obtaining bi- 
partite Ramanujan graphs of arbitrary degrees has appeared recently in |25j . 

Eigenvalues of A on a (q + l)-regular graph may be parameterized by their 
"spectral parameter" s thanks to the relation 

(1.5) A=^-cos(slnq). 

q + 1 



The case s £ K corresponds to A € 



, which is the tempered spectrum. 



9+1 ' 9+1 

In that case we will usually choose s E [0, r] (r = j^t-t)- The case s £ i(— 1/2, 1/2) + 

^Tntl $ e ^) corresponds to A e [— 1, 1] \ I — £r, -^ ) , which is the untempered 
spectrum. The result of this paper will only be of interest in the tempered part of 
the spectrum. 

In what follows, (r n ) will be a sequence satisfying r„ + 2 < R n and q r ' l a n — > 0. 
The sequence (S n ) will be assumed to satisfy S^r^~ 4 — > +oo, for some integer K. 
We also assume that 8 n — > 0, although it is not necessary for the general proof of 
section 0J 

Our aim is to prove the following : 

Theorem 1.3. Let (G„) be a sequence of (q+1) -regular graphs, G n = (V n ,E n ) 
with V n = {1, . . . , n}. Assume that (G n ) satisfies (EIIR) and (EXP). Fix sq € (0,r) 
and let I n = [sq — S n ,so + S n ]. Call (s" , . . . , s™ ) the spectrum of A on G n , and 
(ipi , ■ • ■ , ipn ) o, corresponding orthonormal eigenbasis. 



LetN(I ni G n ) - {je{l,...,n},4 n) e/„} 



be the number of eigenvalues in 



InU Finally, let a n : V n — > C be a sequence of functions such that 



y^ a n (x) = 0, sup|a n (x)| < 1. 

x€V n X 



Note that with the assumptions on 5 n , N(I n ,G n ) — > +00 (see Corollary |5.2l l. 
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The 



Vi £ M n w? 



N(I n ,G n ) 



-> 0. 

>+oo 



Remark 1.4. If a„ does not have zero mean, then by applying the theorem to 
a„-a~, (where ~a~=\ J2xev n a n( x )) w e obtain 



1 V-v I , ,fWl ,ln\, 2 _ 



W(7 „, GJ £ |<tfUtf»>-* 



*-+oo 



Remark 1.5. If we exclude the case so = T /2 in Theorem 1 1.31 we can assume, 
instead of (EXP), the following weaker condition : there exists j3 > such that the 
spectrum of A on L?(G n ) is contained in {1} U [— 1, 1 — 0] for all n. In particular, 
the theorem applies for bipartite regular graphs in this case. 

We can also say something in the case sq = t/2 for bipartite expander graphs, 
that is if there exists /3 > such that the spectrum of A on L 2 (G n ) is contained 
in{— 1, 1} U [— 1 + (3,1 — j3\ for all n. We need to strengthen the condition on the 
functions a n (x) in the theorem for the conclusion to apply : if V n — V„ U V 2 is the 
bi-partition of V n , then we need that 

^ a„(x) = ^ a n( x ) = °- 

The theorem then tells us that we have equidistribution of most eigenfunctions with 
eigenvalue near t/2 on each set \Q and V„, without providing information on the 
relative weight of these two sets. 

Remark 1.6. The proof will show that we can weaken condition (EXP), by 
allowing the spectral gap j3 to decay with n "not too fast" (/3 3> r„ ' is enough). 
See Section |H1 

Since random (g+l)-regular graphs satisfy both (EXP) [29lE[IH and (EIIR) 
[27) , our theorem applies to them with the values of R n , a n given in Example [T] 

Corollary 1. Let a n : V n = {1, . . . , n) — > C be a sequence of functions such 
that 

} J a„(x) = 0, sup|a„(a;)| < 1. 

l£V„ X 

Choose (G„) uniformly at random amongst the (q + l)-regular graphs G n — 
(V n , E n ) such that V n — {1, . . . , n}. Choose j uniformly at random in N(I n , G n ). 
Then for any hxcd e > 0, 



P 
when n — > +oo. 



( (^™ ) ,a„V. 



3 I 



>e — ^0 



The statement of Theorem ll.3l is the exact analogue of the Shnirelman theorem 
in its form (|l.lj) . However, we do not have a statement analogous to the convergence 
of measures (|1.2p , because our sequence of measures does not live on a single space; 
instead, it is defined on the sequence of graphs G n . We do not know of a notion 
that would be adapted to describe the limit of the family (G n ) endowed with the 
probability measure (\tp) n (x)\ 2 ) xe v„. 
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We can generalize Theorem ll.3l bv replacing the function a with any finite range 
operator : 

Theorem 1.7. Let (G n ) be a sequence of (q+1) -regular graphs, G n — (V n ,E n ) 
with V n = {1, . . . , n}. Assume that (G„) satisfies (EIIR) and (EXP). Fix sq £ (0, r) 
and let I n = [sq — S n , sq + 5 n ]. Call (sj , . . . , s„ ) the spectrum of the laplacian on 
G n , and (ip\ n , . . . , ipn ) a corresponding orthonormal eigenbasis. 

Let N(I n , G n ) = \{j € {1, . . . , n}, sj° £ /„} . 

Fix DgN. Let A n be a sequence of operators on L 2 (G n ) whose kernels K n : 
V n x V n — > C are such that K n {x, y) = for d(x, y) > D. 
Assume that sup x yf zy n \K n (x,y)\ < 1. 
Then there exists a number A n (so) such that 






i 

— > 0. 

n )-+oo 



With the notation of ^S. 6 2[ we can write A n — Op(a n ), a n £ S® ; and we have 
the expression A n (s ) = i J2x£D n In a n{x,U, s )dv x (u>). 

Quantitative statements (i.e. rates of convergence) will be given in Section [6) 



2. Theorem 11.31 : outline of the proof in the case so = r/2 

We first give a proof in the special case So = T /2- The reason for treating this 
case separately is that one can give a proof which is exactly parallel to that of the 
Shnirelman theorem on manifolds [3i2j [6l 1361 137] . The case of arbitrary Sq requires 
additional arguments and will be treated in Section |4j 



2.1. Upper bound on the variance. Fix an integer T > 0. Let x be a 

smooth cut-off function supported in [—1,1] and taking the constant value 1 on 

[-1/2,1/2]. We write 

(2.1) X.(.)-x(^? 

so that Xn = 1 °n I n - We use the pseudodifferential calculus and the notation 
defined in Section |3l taking the cut-off parameter r equal to r n (from condition 
(EIIR), as explained before the statement of Theorem ll.3l) . 

To simplify the notation, we will write ipj = %/)j, Sj = s" , and a = a n . 
The observable a is a function on G n , in other words a T„-invariant function on 
X. Let il be the boundary of X (see section [3J , then a extends to a function on 
X x il x [0, t] that does not depend on the last two coordinates. The notation Op(a) 
is then defined in section [3j 
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Thanks to Lemma 13.101 and to the "Egorov property" Corollary I3.9| we havq3 

n n , _. 

X]x™( s i) 2 |^' a ^)| 2 = ^IW^OpcJaXrOV^^ + ^O ( — — 
3=0 3=0 ^ [ n) 

+0(T 2 q r+2 )\{x G D n , p(x) <r + 2}\J j£(a)d/i(a) 
+0(T 2 )n I 'o(\s - s \ 2 )xl(s)dn(s) 



where a T := i J2k=o a ° °" 2fc an( ^ c : X x £1 — ► X x fHs the shift (see 
Next, we use Lemma 13.31 to write 



^2\{i>j,OPG n ( aT Xn)i>j)\ < \\Op Gn (a T Xn)\\ 2 Hs 

3=0 



- 5Z / W T (x : uj)\ 2 diy x (uj)xn(s) 2 dp(s) 

+q r ^2 / \a T (x,uj)\ 2 dv x (uj)xn{s) 2 d[i(s) 

x£D rl ,p(x)<~r 

< J2 f \a T {x,uj)\ 2 diy x (uj) [ Xn(s) 2 dn(s) 

x£D n 

+q r \{x G D n> p(s) < r}| / Xn(sfdn(s). 



We also know from the Kesten-McKay law (Section [5j Corollarv l5.2p that 



X n(s) 2 dp( S )=0(l). 



N{I n ,G r , 



5 To prove the extended Theorem 1 1,71 we also need Lemma l3.11l 
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We thus have 

— J-— £ \(^,a^)\ 2 = 0(l)- J2 f\a T (x,w)\ 2 du x (u;) 

+ 0(T 2 q 



2^r+2M X ^ D n,p{x) <r + 2}| 



" ^0(J—) + 0(^) +( >iT^ i} 



N(I n ,G n ) \(r8 n )°°J \r 

0(1)- £ [\a T (x,u)\ 2 dv x (uj) 

+ 0(T 2 q r )a n + 0(1)^-0 ' 



<5» \K) C 
T 2 / 

2 



O ( ^ ) + 0(2*£) 



Our choices of r = r n and 5 n imply that the last four terms vanish as n goes to 
infinity while T is fixed. 



2.2. Expansion and ergodicity. We write, using 

T— 1 T— 1 

For every ifl and fceff, define the partition of fl in (q + l)q k ~ 1 setqj 

f2(x,y) := {w S fl | [»,w) = (#,#1,0:2, ■ ■ ■) and a; fc = y} , 

where y € X and d(x,y) — k. Then w i— > a o <r (#, w) is constant on f2(rc, y) for 
every y G X such that d(x,y) = k, and 1/3. (fi (a;, j/)) = , +1 j ^_ t . We have 

\J / aod'fijwjafijdi/jfw) = \J 2_. / ao cr fc (x,aj)a(x)(iz/ 2 ;(a;) 

dx{x,y)=k 

= 22 S k a(x)a(x), 
xeD n 
where So — Id and for all k > 1, S k is the stochastic operator defined as follows by 
its kernel on the tree X : 

( 2 - 2 ) S k f(x) = - ] fc _ 1 £ /(»)■ 

We thus have 

T— 1 T— 1 

~ E / |a T (a:,w)| 2 d^(w) = -^ ^!E E S 2 |k-j|a(#)a(a;). 

On the quotient G n , the spectrum of S k is the set {4>sj (k),j = 1, . . . , n) where 
4> s is the spherical function, 

(2.3) *.(*) = ^ (4- cosiksHq)) + i_l ^ + |Mn( g ) A 

V 9 + 1 9+1 sin s ln(gj / 

See section [3] for a definition of the boundary Q and of the notation [x,u>). 
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and {sj,j — 1, . . . ,n} are the spectral parameters for the operator A defined by 

dm. 

Using the parameterization (|1.5[l of the spectrum, the eigenvalue A = 1 corre- 
sponds to 

1 

is = — . 
2 

Because of the (EXP) condition, the other untempered eigenvalues satisfy is £ 
(0, 5 — j3) or is + ijj- £ (0,^-/3), for some /3 > independent of n. It follows that 
\4>s -{k)\ < Cq~P k with C,(3 independent of n. The eigenvalues of the self-adjoint 
stochastic operator ^ X)fc=o Ej=o &2\k-j\ are therefore bounded by 

T-1T-1 
_ V V „— 2/3|fc— _7| < ° 

T 2Z^Z^ y "2T(l-q- 2 ' 3 ) 



in modulus. They are contained in {1} U ~Yg) Sa for some C, independent of 
n, T, (3 (the eigenvalue 1 has multiplicity 1, corresponding to the constant function). 
Thus, if a satisfies YlxeV a ( x ) = an< ^ su Px \ a i x )\ — Is we have 

1 1 T-1T-1 



/ 1 \ 

7, — n -; — n ™./- n \ r* / 



nT 2 

k=0 j=0 xG-D 



2.3. Conclusion. We obtain, using the results of the previous sections and 
the Kesten-McKay law fCorollarv l5.2|) . 



' £ \^M 3 )\ 2 <0(— ( -J— )+o(^)+0(THl) 



N(I n ,G n ) trl " " Jn ~ \S n {rS n )-/ \r 

+ 0(T 2 q r a n ) + 0(^ 

3 

If we choose the sequences r = r n and S n satisfying q r a n — > and s Z s rg — > 
for some integer K, we finally have 

As the left-hand side of the equality does not depend on T, we take the limit 
T — > oo to obtain 

n^oo N(I n ,& n ) ^ 

3. Elements of pseudodifferential calculus 

In JJ3.1H3.2I we recall some of the tools of pseudodifferential calculus that were 
introduced in [23] . However, the following important remark has to be made : in 
order for Theorems 11.31 and 11.71 to have full strength, we should not impose on the 
symbols a too strong regularity conditions, that would have the effect of making 
the theorems trivial consequences of the (EXP) condition. Thus, we pay attention 
to only use from 1231 the properties that do not require regularity of a with respect 
to the x-variable. 
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In the following sections we try to construct a pseudodiffercntial calculus on 
the quotient. 

3.1. Definition of Op(a) on the infinite (<7 + l)-regular tree. Let il be the 

boundary of the tree. It is the set of equivalence classes of infinite half-geodesies 
of X for the relation : two half-geodesies {xi , x% , X3 , . . .) and (2/1 5 J/2 ? S/3 5 •• •) are 
equivalent iff there exist k,N £ N such that for all n > N, x n +k = Vn- For every 
u) £ fi, we will denote by [x, u>) the unique half-geodesic starting at x and equivalent 

tO L). 

Let a : 1 x O x [0, t] — > C be a bounded measurable function. In |23j . the 
operator Op^(a) was defined by 

Op x (a)u(x) = V f f q^ +ls ) {h - {v) - h - {x)) a{x 1 Lo 1 s)u{y)dv x {u)d i x{s) 
yeX JnJo 

for every u : X — > C with finite support. Here h u (x) is the height function (or 
Busemann function), d[j,(s) is the Plancherel measure associated to the (q + 1)- 
regular treqj, and v x is the harmonic measure on 51, seen from the point x. We 
refer to [7] for more background. We will denote by 

K x (x,y;a) = f /V**') (Mw) ~ Mx M*,«,«)^(w)^(*) 

Jn Jo 

the kernel of Op x (a) . 

3.2. Class of symbols. From [7], we know that the fact that Kx(x, y;a) = 
for dx(x 7 y) > D is equivalent to the four following conditions on a : 

• a is a continuous function; 

• a extends to a 2r-periodic entire function of exponential type D uniformly 
in w; i.e. for all x there exists C(x) > such that 

\a(x,w,z)\ < C(V)g D|am(z)l Vw G n, Vx e X; 

• a satisfies the symmetry condition 

q^- is ^ h ^^' h ^ x ^a{x,0J,s)dv x {u) - / q^ +ls){h - {v) - h -' {x)) a{x,u J ,-s)dv x {u J ) 
m Jn 

for all s £ C, x e X. 

• a is a D-cylindrical function, that is : if the two half-geodesies [x,ui) = 
[xq,Xx,X2, ■ ■ ■) and [x 1 ,lo') = (x'q, x' x , x' 2 i ■ ■ ■) satisfy Xj = x' for < j < 
D, then a(x, u, s) = a(x' , u>', s). 

We shall denote by S®(X) the class of such functions. In [23] . another class of 
symbols was considered : 

Definition 3.1. S(X) is the class of functions a : X x fl x [0, r] — ► C such 
that 

• for every (x, w, s), for every k G N, d^a(x, w, s) is continuous on X x 51 x 
[0, t], and for every Z G N, there exists C; > such that, for all n £ N, for 
every (x,o;,s), 



|(o-£^o)(s,w,*)| < 



(1 + n) 1 



'dfjb(s) = \c(s)\ 2 ds where c is the Harish-Chandra function. 
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• for every k € N, and every (x, w), d^a(x, w, 0) = c^ a(s, w, r) = 0. 

Here £%a(x, u>, s) is the projection of a(x, w, s) on functions depending only on 
the first n vertices of the half-geodesic [x,ui) (see |23j for a formula). In particular, 
£*a(x, s) — a(x, s) if a does not depend on u>. 

It is proven in [23 that S(X) endowed with usual addition and multiplication 
is an algebra. This makes it more suitable for semiclassical analysis than the class 
S„(X). It also has the property, crucial for us, that 

a £ S(X) =>■ a o a G S(X) 

where a is the shift, a(x,u) — {x\,ui) if [x,uj) = (xq,x\, X2, ■ ■ •)■ It is proven in 
[23] that 

/ M+l \ d(x,y) 

(3.1) \K x (x,y;a)\<C[\\a\\n,M+Yj\d*a\U\ ' 



for any M, where ||a||n,M = sup(x,u),«) su Pn(l + n) M \a — £*a\(x,uj, s), and that as 
a consequence, Op x (a) extends to a bounded operator on L 2 (X) if a € S(X). 

If a(x, oj, s) — a(x) depends only on x, then Op x (a) is the operator of multi- 
plication by a. At several places we will use the fact that 

(3.2) Op 3E (o)Op 3E ( V )=Op 3E M 

if (p = y(s) only depends on the last variable and a(x,a;, s) € S^X), say. 

In most of what follows, we will actually need very few conditions on the sym- 
bols a(x, uj, s). Essentially it will be required that a € L°° and that for every x £ X, 
a(x, •, •) S L 2 (f2 x [0, r]). For convenience, we can assume that a € Sg(X) for some 
D € N, or that a € S(X). In Lemma [3.111 we use the condition a G S„(X). 

3.3. Definition of Op G (a) on a finite graph. Recall that G„ is written as 
a quotient T n \X, where T n is a group of automorphisms of X, whose elements act 
without fixed points. 

Let us now assume that a is r n -invariant, meaning that a(j • x, 7 • u), s) — 
a(x,to,s) for all (x,ui,s) and all 7 £ T n (where the action of T n on the boundary 
Q is obtained by extending its action on X). For a r n -invariant symbol, we have 

K x (~f-x,~fy;a) = K x (x,y;a) 

for all x,y € X and 7 € T„. The proof of this fact is identical to the proof of 
Proposition 1.1 in [35 . 

We now define Op(a) on the quotient. 

Definition 3.2. Assume the sequence (G„) satisfies (EIIR). Let r = r n be a 
positive number. 

If a is r„-invariant, we define Op G (a) to be the operator with r n -bi-invariant 

kernel 

TJ . . . v-^ r , . . / d(x,j ■ y) 

K Gn {x,y] a )= 2^ K x(x,i -y^jx I 

7er„ ^ 

Here x is a cut-off function that satisfies the conditions of ij2.ll (although it 
need not be the same cut-off as in £|2. 1 1 we use the same notation). 

Compared to the case of manifolds, a difficulty we meet is that we are not 
able to prove that Op Gn (a) is bounded on L 2 (V n ) independently of n (actually, 
inspection of simple examples show that our conditions on a are not sufficient 
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to ensure this). Note however that we are only interested in Op Gn (a)ipj for 
Ay 1 in the tempered spectrum; more precisely, we shall only need to estimate 



quantities such as 



N(I n ,G„) 



£., 



'ei„ 



dn) 



■An), 



For that purpose it 



will be sufficient to know that the Hilbcrt-Schmidt norm of Op G (a) does not grow 
too fast : 



Lemma 3.3. We have for every r > 

|0 Pg »||^< J2 E \Kx(x,y;a)\ 2 - 
xeD n y ex 

p(x)>r d(y,x)<r 

- E / \a{x,co,s)\ 2 du x (uj)dn(s) + q r 

l6D„ 

Proof. By definition we have 

\\Op Gn (a)\\H3 = E \ K G n {x,y;a)\ 

x,y£D n 



E 

xeD n v ex 

P(x)<r d(y,x)<r 



E 

x€D n 

p(x)<r 



]T \K x (x,y;a)\ 2 
a(x,u>, s)\ dv x (u)dp,{s) 



E 

x,y£D n 



^2 K x(x,J-y;a)x 

76T„ 



d{x,i-y) 



We split the sum into two parts, whether p(x) > r or not. If p(x) > r, then the 
sum over 7 G T„ is reduced to only one term, thanks to the cut-off function. If 
p(x) < r, then there are at most q r terms in the sum over 7 e r„, and we can use 
Cauchy-Schwarz inequality to bound it as follows 



|0 PG »||^ ^ 



< 2, l_ \K X (x iy ;a)\ 2 + q r J2 E l^e(x,y;a)| 5 

xeD n y£X xGD n y££ 

p(x)>r d(y,x)<r p(x)<r d(y,x)<r 

< E Ei^.»;«)i 2 +« r E J2\ K *( x >v> a )\ 2 - 



x£D n yGX 
p(x)>r 



xeD n yGX 
p(x)<r 



Plancherel formula for the Fourier- Helgason transform, applied to y 1-4 K%(x, y; a) 
with a; fixed, converts this last expression to 



\Op Gn {a)\\ 2 HS < ^2 Hx,uj,s)\ 2 dis x (uj)dp(s) 

|a(x,cj,s)| dv x (uj)dp{s) . 



16A 



E 

p(a;)<r 



a 



3.4. "Egorov"-type properties. For Quantum Ergodicity on manifolds, the 
Egorov theorem is a statement saying that the matrix elements (ipj,Op(a)tl>j) re- 
main almost invariant when transporting a along the geodesic flow, when ipj are the 
eigenfunctions of the Laplace-Beltrami operator A. This is proven by showing that 
taking the bracket [A, Op(a)] amounts to differentiating a along the geodesic flow 
(up to some "negligible" error term). Here we try to perform a similar calculation. 
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We define a map a : X x fl — > X x £1 by a(x,u>) = (x',ui) where x' € X 
is the unique point such that dx(x,x') — 1 and x' belongs to the half-geodesic 
[x,u>). If [x,lo) — (x,X\,X2,X3, . . .), then a(x,u>) — (x\,u>), corresponding to the 
half-geodesic [x\,w) — {x\,X2,Xs, ■ ■ •)■ I n symbolic dynamics, a is the shift, and 
if we compare with Quantum Ergodicity on manifolds, a plays the role of the 
"geodesic flow" on phase space IxO. This map is not invertible, actually each 
point has exactly q pre-images. We shall denote by U the operator a \— > a o a. For 
a : X x il — > C, we define La : X x O — > C by 



La(x,ui) = - ^2 a (v^). 



q 
If a and b are compactly supported functions, we have 

y / a o <t(x, uj)b(x, Lo)dv x (m) — > / a(x, u))Lb(x, uj)dv x (w), 

in other words L is the adjoint of U on the Hilbert space L 2 (X x $7, ^^ ^^^(w)). 
In addition, we also have LU = I, reflecting the fact that U is an isometry of L 2 (Xx 
0, ~^2 x 5 x dv x {uj)). The operators U and i preserve the T„-invariant functions. If a 
and 6 are a r„-invariant functions, we still have 

(3.3) \ / a o <r(x,uj)b(x,uj)di' x (uj) = N / a(x,ui)Lb(x,uj)di , x (uj). 

xG.D„ ° xeD n , ' n 

Recall that D„ is a fundamental domain for the action of T n on X. 

In what follows, we extend the definition of U and L to functions on 1 x O x K, 
by a trivial action on the last component. The crucial bracket calculation is the 
following : 

Proposition 1. If a e S(X) is r„-invariant, then 

[A,0p Gn (a)] = 0p Gn (c) + R 

where c £ S(X) is a symbol given by 

q 1/2 
c(x, U>, s) = (q ts (a o a — a)(x, ui, s) + q ls (La — a)(x, ui, s)) , 

and R is a remainder such that 

\R\\ 2 hs<o(^J I J2 j \<x,u,s)\ 2 dv x {w)dn(s) 

^ / \a(x,uj,s)\ 2 dv x (uj)dn(s) \. 



.x£D n 

„r+2 



Q 



xGD n 

p(x)<r+2 



Since {ipj, [A, Op G {a)\ipj) = for every laplacian eigenfunction ipj, this implies 
(3.4) ^|(Vi,Op Gn (c)^)| 2 <0 (4) ( E f\a(x,uj,s)\ 2 du x (uj)dn(s) 

^ / |a(a;,cj,s)| 2 d^ :E (tj)d^(s) . 



„r+2 



i6D„ 
p(:r)<r+2 
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Remark 3.4. An exact analogue of the usual Egorov theorem on manifolds 
would require an estimate of ||Op(c)|| ^,2 and show that it is up to a vanishing 
remainder term. Here, due to our use of the Hilbert-Schmidt norm, we will only 
show this for the average 

sf'ei„ 
which is sufficient to prove our theorem. 

PROOF. Let us denote by K Gri (x,y; [.]) the kernel of [A, Op Gn (a)]. We know 
from |23j that Kx(x,y;c) is the kernel of [A, Op x (a)]. We are interested in the 
difference K Gn (a;, y; [.]) — Kg„ (x, y; c), which is the kernel of the operator R. 

We have 

K G„(x,y,[.\) = — -I Yl K G n ( z ,y,a)- X! K G n (x,z;a)\ 

5 \d(x,z) = l d(z,y) = l J 

(3-5) =r^y E^s(a:,7-»), 



q ■' -7^ 



where 



(3.6)^(x,y):= ]T ff s (*, y; a) X (^- ) - £ #*(*,*; a)* (^^) 

d{x,z) = l ^ ' d(z,y) = l ^ ' 

Because of the cut-off functions, the sum (|3.5j) only runs on those 7 € r„ for which 
d(x,7- y) < r + 1; and in (f3~6]> we have if x (a;, y) = J^3e(«,l/)l{d( a ;, v )<r+i}- 

In the first sum of the right-hand side of equality (|3.6|) . d(z,y) = d{x,y) ± 1, 
because x and z are neighbours. In the second sum d(x, z) — d(x, y) ± 1, because z 

and y are neighbours. Since \ is a smooth function, both x f ^ v ' ) and % I ^ x ' z ' J 
are equal to \ I r ) "*" ^ ( r ) » anc ^ we nave 



^(x,y)=( J^ Kx(z,y;a)- J^ ^(a;,z;a) J x 

iii(x,z)=l d(z,y) — 1 



d{x,y)' 



° ( 7 ) ( X! I^(^,y;a)|+ X! 1-^(3, 2; a)| J fl{d(x,2/)<r+i} 
'd(x,y)~ 



d(x,z) = 1 d(z,y) = l 



(q+ l)K x {x,y;c)x 



Of-) J^ 1^36(^,2/; a)| + X! I^( a; > z ; a )l ]1 {d(x^)<r+l}- 



T , , 

d(x,z) = l d(z,y) = l 
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Now if we go back to (|3.5|) we get KG n (x,y; [.]) = KG n {x,y;c) + KG n (x,y;R), 
where Kc n {x, y; R) is the kernel of the operator R, given by 

/ \ 



K Gn (x,y;R) = 0(l) E E \Kx{z^-ya)\ + E \ K ^ a 



7er„ \ zex 

\d(x,z) = l 



zex 

d(z,i-y) = l 



1 



{d(x,f-y)<r+l}- 



We estimate the Hilbert-Schmidt norm of R by first writing 



E E 

i£D„ yG-D„ 



E E 



\ 



E E \ K x(z,^-y,a)\ + E \K x {x,z;a) 

7Gr„ I 2 £l z£X 

\d x (x,z) = l d x (z,j-y) = l 



11 



{d s (a;,7-i/)<r+l} 



/ 



E E l- ftr *'( z '7-2/;a)l ] l{d^(z, 7 -y)<r+2} + E \ K x(x,z;a)\i {dxix ^. y) < r+1} 



76T„ z£j 



d x (z,j-y) = l 



We then use Cauchy-Schwarz and reason along the same lines as in lemma [5751 to 
bound the former expression by 

<2(« + l)( E E E \Kx(z,y,a)\ 2 + E E E l**(a:,*;a)| 2 

V z£-D„ y£l a;e_D„ xED n y£X z<£X 

p(z)>r+2 da n (z,x)=l p(x)>r+l d x (z,y) = l 

+ f +2 E E E \Kx(z,y,a)\ 2 + q r+1 E E E l^e(w)| 2 ) 

zeD n y£X x<ED n xeD n y£X z£X / 

p(z)<r+2 d Gn (z,x) = l p(x)<r+l d x {z,y) = l 

<2( q + i) 2 ( E £!**(*» w«)l 2 + E £l#s(s,*;a)l a 

V 2£fl„ y£X x£D n zGX 

p(z)>r+2 p(x)>r+l 

+i r+2 E Y,\ K ^ z ^ a )\ 2+( i r+l E £!**(*, *;a)i 2 ) 

z£fl„ y£l i6D„ z£l / 



p(z)<r+2 



p(x)<r+l 



< 



4( 9 + l) 2 E \a{x,uj,s)\ 2 dv x (Lo)dn(s)+q r+2 E j \a{x,u,s)\ 2 dv x (uj)d^(s). 



i6D„ 

Finally we have 



xGD n 
p(x)<r+2 



\R\\ 2 hs<o^ ( E j Hx^,s)\ 2 d Vx {ujW{ s ) 



,i£D„ 



a;G-D„ 
p(x)<r+2 



D 
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In what follows, Proposition [1] will be translated into an invariance property of 
the type 

A^k) 5: K^,Op G „(a)^.)| 2 ~ j^^ g E l^,Op G „(Ta)^>| 2 

for some operator T. A key idea is then to take advantage of the spectral properties 
of T and its iterates T k for k <G N. In the special case where so = t/2 we can take 
T = U 2 fCorollarv l3.9p , which makes this spectral value special. For general values 
of so, T is a linear combination with complex coefficients of the non-commuting 
operators L and U, and its spectral properties are not so nice. The aim of the 
successive operations done in Corollaries 13.51 to 13.81 is to replace Ta with Ua up to 
some error term. We first replace a with q ls a o a in (J3.4I) to obtain 



Corollary 3.5. 



<',->l 2 



Y^m,Op((U-I)(q 2is U-I)a)^ 

v 7 i£D„ J 

+ 0(q r )\{x G D n ,p(x) <r + 2}\ [\\a(; •, s^d^s) 



- ° ( 62 ) ^ / \a(x,uj,s)\ 2 dv x (uj)diJ,(s) 
xeD n 



where Ua = a o a 

Proof. Recall that symbol c of proposition [1] is given by 

„l/2 

' (q ts (Ua-a)+q- ls (La-a)). 



q + i 

If we replace a with q ls Ua we have 

c = ^ (g 2 "'(t/ 2 a - Ua) + {a - Ua)) = S. (TJ - I)(q 2 ™U - I)a. 

q+ 1 q+ 1 

It follows from (157311 that 

J2m,Op((U-I)(q 2is U-I)a)^)\ 2 
3 

<o(— J ^ \Ua(x,uj,s)\ 2 dv x (uj)dfi(s) 

\ r ' x£D n J 

+ ($) J2 f \Ua{x^,s)\ 2 dv x {uj)d^{s) 

<0(— J E \a(x,uj,s)\ 2 diy x (uj)diJ,(s) 

+ ° (^) K* e D ». rts) ^ r + 2 >l /lM. •> s )\LMs), 

where we used the fact that U preserves the L 2 and L°° norms. □ 

The idea is then to invert (q 2ls U — I). As the series J^k q 2%ks U k is a formal 
inverse to (q 2ls U — I), we apply Corollary 13.51 to a = X^^o q 2lks U k b := &jv-i, 
where 6 € S'(X) and AT is an arbitrary integer. We obtain 
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Corollary 3.6. 

Y,\{^^{{U-I){I-q 2iNs U N )b)^)\ 2 

3 

+ O (¥f) \{x € D n , p(x) <r + 2}\ [\\b(; -, S )||^dM(s) 



PROOF. We apply Corollary [331 to a = ^aLc/ q 2lks U k b := &jv-i and use the 
identity 

(q 2ls U - I)b N -! = (q 2iNs U N - I)b 
combined with the fact that [/ preserves the L 2 and X°° norms. □ 

If we apply Corollary 13.51 to a = -k X}fc=o ^*j we obtain 
Corollary 3.7. 

Y,\(^, Op (ub-b- q 2 ™u(u-i)b^)^)\ 2 

^2 / |6(^,w,s)| 2 d^(a;)(i/x(s) 



,'iV 2 
<0 



KG-Dn 



i ^~f ) K* e ^n.pCs) < r + 2}| / ||&(-,-, S )||^( S ) 



Note that the "remainder term" q 2ls U(U — /)fo( s,Ar ) is not small in the symbol 
norm : in Section [H the (EXP) assumption will be used to show that it is small 
in the L 2 -norm. This is a major difference with the Egorov theorem on manifolds, 
where no ergodicity assumption is needed. 

Proof. We know from the proof of the previous corollary that 

(I - q 2ts U)b k = (I - q 2 ^ k+1 ^U k+1 )b. 

It follows that (I - q 2is U)± J2k=o bk = b- q 2i8 Ub^ 8 ' N \ and 



N 

N-X 

N) 

N 



N-l 

(u - /)(/ - q 2ts u)— J2 b k = Ub-b- q 2ls U(U - I)b^ s 



k=Q 



□ 



Combining with the Hubert-Schmidt estimate, Lemma l3.3[ we get 
Corollary 3.8. 



Y^K^Pj, Op (Ub - b) ^)\ 2 <o(^\ J2 l\b{x,w,s)\ 2 dv x {u)dix{ S ) 

j \ r / x £D n •* 

+ 2J2 [\b (s ' N) (x,u},s)\ 2 dp x (Lj)d»(s) 



x€D n ■ 

+ 0(N 2 q r )\{x G D n ,p(x) < r + 2}| / ||6(-, ■, s^d^s) 
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PROOF. We write 



J2 Mi, Op (Ub - b) ^-)| 2 < 2 J2 I (1>j, OP (^ " & " 9 2 '^(^ - /)& (s,JV) ) il>j) 

j j 

+ 2^1 (Vj, Op (q 2ls U(I - U)b^ s ^ N A ipj 



The first term on the right-hand side is estimated by Corollary 13.71 We estimate 
the last term thanks to Lemma [3. 31 and we use the fact that U preserves the L 2 
norm by U : 

J2\{^,Op {q 2is U(I-U)b^^)\ 2 

3 

< E I \(I -U)b^ N \x,u,s)\ 2 dv x {uW{s) 

+ q r J2 J\(I-U)b^ N Hx,u;, S )\ 2 du x (u;)diJi(s) 

x£D n ,p(x)<r 

< J2 f\b^ N \x,oj,s)\ 2 du x (uj)dn(s) 



xeD„ 



+ AN 2 q r \{x e D n , p(x) < r}\ / ||6(-, •, s)\\i,d^ 3 ) 



a 

As already mentioned, the value so = r/2 is special and the previous corollaries 
may be replaced by the following, simpler one. In case the support of a shrinks 
around So = r/2, this is a closer analogue of the Egorov theorem on manifolds in 
the sense that no ergodicity or expanding assumption is needed to show that the 
remainder term goes to 0. 

Corollary 3.9. We have 
^K^OpM-OpCaoa 2 )^)! 2 ^^ Y, f\a(x,uj,s)\ 2 diy x (cj)dn(s) 

+ O (q r+2 ) \{x € D n , p(x) <r + 2}\ J\\a(; •, s)\\ldn(s) 

+ n fo(\ S - Sa \ 2 )\\a(;; S )\\ldt,( S ) 



Proof. We replace the symbol a in proposition[T]with q ls aoa. As L(aoa) = a, 
the symbol b becomes b = a — Ta, where 



Ta(x,u, s) = —q 2ts a o <t 2 (x,oj, s) + 2q ts cos(slng)a o a(x,u!,s) 
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and we have J2j KV'j; Op(a) — Op(Ta)tpj)\ 2 < II^IIhs> where 
\\R\\ 2 HS < O (^) ]T f\aoa( Xl uj, S )\ 2 d Vx {oj)d^ S ) 

+ 0(q r+1 ) ^ \aoa(x,u,s)\ 2 dis x (uj)dii(s) 

x£D n ,p(x)<r+l 

<Oi— J ^2 \a(x,uj 7 s)\ 2 dv x (uj)dfi(s) 

\ T ' x£D n •* 

+ O (q r+1 ) \{x G D n , P (x) <r + l}\ f\\a(; -, *)||^d/i(*)- 

Now Ta — a o a 2 — c, with 

c(x,a;,s) = (1 + q 2%s )a o er 2 (x,w,s) + 2g' s cos(slng)ao a{x,uj,s) 
and wc can write 



^|(^,Op(a)-Op(aoa 2 )^)| 2 <2^|(^,Op(c)^)| 2 +2||i?! 

3 3 

<2\\Op{c)\\ 2 HS +2\\R\\ 2 HS . 



2 
HS 



Recalling that ,sn = „? , wc have 

° u 2 In g ' 



|Op(c)||#<j< Yl / K 1 + « 2lS ) a ocj2 + 2 1 lS cos(s\nq)aoa)\ 2 dv x (uj)dfi(s) 
xeD n •* 

+ q r Y \(1 + q 2is )aoa 2 + 2q is cos(s\nq)aoa)\ 2 dv x (u})dLi(s) 

x£D n ,p(x)<r 

< 2_, I 0{\s — so| 2 )(|a o a 2 \ + \a o a\) 2 {x 1 cu, s)dv x (uj)dii{s) 
xeD„ J 

+ q r 2J 0(\s — s a \ 2 )(\ao a 2 \ + \ao a\) 2 (x,LU,s)dv x (Lu)d/j,(s) 



xGD n ,p(x)<r * 

<n I 0(\s — s \ 2 ) sup \a(x,ui, s)\ 2 d/i(s) 



+ q r \{x E D n ,p(x) <r}\ / 0(\s - s a \ ) sup \a(x,u, s)\ dfi(s) 

X,U) 



a 



3.5. Two more formulas about Op G (Xn)- 

Lemma 3.10. 



with 



OpaMn)^^^^ 



\f ] =Xn( Sj ) + r*0< ' 



if Sj € [0, r] (tempered eigen] 'unctions). 
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in) 

Proof. First note that ijy- is associated to a r„-invariant eigenfunction of 
the laplacian on the tree X, that we will still denote by ipj 1 ■ We have on the tree 

Op Gn (xn)4 n) (x) = £>*(*,!,; x»)x (^y£) ^\v) 

yex \ r / 

and Kx(x, y; Xn)x ( r ' v ) depends only on d(x, y) because \n does not depend on 
(x,uj). We thus have 

where f(sj) is given by the spherical transform of the kernel 

f( Sj ) = ^K x (x,y; X n)x (^^-) <t> Sj (d(x,y)) 

vex \ r / 

and <p s . is the spherical function associated to Sj defined in (|2 .3[) . Now 

f( a i) = ^2 K x(x,y;Xn)<p S] (d{x,y)) 
vex 

'd(x,yY 



- Yl K x(x, y; Xn)(t> Sj (d(x, y))(l-x 
vex ^ 

Xn(Sj) 

J2 K x (x,y ]X n)ct> s Ad(x,y)) (l- X ^' U ' " ] 



vex 

Because Xn € S(X), according to the rapid decay property of the kernel of 
pseudodifferential operators (13.11) . we have 

CM+1 \ d(*,v) 

gfcn|UJ (1+ % ( ^ ))M 

for every M S N. Moreover, if Sj is a tempered eigenvalue, then 

So, using also the fact that ||f?gXn||oo < CkS^ k , we obtain that for every MeN 

\f(s)-x MK^y q ~ d(X ' y) (i-x(^ yY 

y t x. 
Cm+1 V q -d(x,y) 

- 6™ +1 (l + r) M - 2 ^(l + d(x >y )) 2 
Cm+1 \~y v-^ 1~ 



Cm+i 

M+l(l +r )«»-2 £j (1 + k) 



S S 7T 

y— 



5^ +1 (l + r) M - 2 ^ ^ (l + k) 2 

n ^ ^ I keNy:d(x,y)=k V ' 



We thus have 



l/(ai)-Xn(«i)|=r d O l (r(5n)M+ i 
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for any M . □ 

Lemma 3.11. Fix an integer D. Let a be such that K x (x,y;a) = for 
d x (x,y) > D (in other words, a G S®(X)) and ip = ip(s). We have 

°Pg„ ( a f) = °Pg„ ( a ) °Pg„ (f) + R 
where 

\\R\\ 2 hs < ^(r- 1 ) J ^ 2 (s)d/x(s) (n + \{x € V n , p(x) < r + D}\q r+D ) 
PROOF. The kernel Kg„ (x, z\ aip) of Op G (aip) is obtained by the periodization 



d(x,7 • z) 



K Gn (x,z;aip) = ^ #a:(a;,7 • 2; atp)x 

7er„ 

Because <p only depends on s, we note that 

K x (x,z;aip) = J^ K x (x,y;a)K x (y, z;<p) 
yex 

(see ([3T2])) and 

K x (x, z;aip)x f ! — J = 51 K x (x,y; a)K x (y, z; p)\ f ! — 

^ r ' yea; ^ r 



^^, P M„)(xf^-x /%2: 



= X! ^( x ' 2/' a ) K x(y, z; v)x f — ! — 

yes ^ 

ye* 

After L„-periodization, we note that I] 7 er„ 2»e£ ^(^ 2/; a)K x (y, •y-z; <p)x ( V 'J' Z ' ) 
is the kernel of Op G (a) Op G (ip) (as soon as D < r). Using Cauchy-Schwarz and 
the fact that K x (x, y; a) is supported near the diagonal, the Hilbert-Schmidt norm 
of the operator with kernel 

^2 ^2 \K x (x,y;a)\\K x (y,-f ■ z;tp)\i dx(x ^. z) < r+D 

76r„yG3£ 

on L 2 (G n ) can be bounded by 

q 2D sup \K(x, y; a)\ 2 f v 2 (sW(s) (n + |{a; e V niP (x) < r + D}\q r+D ) . 

x,y J 

D 

4. The proof for arbitrary so 

4.1. Upper bound on the variance. Fix an integer T > 0. Let x be a 

smooth cut-off function supported in [—1,1] and taking the constant value 1 on 

[-1/2,1/2]. We write 

' s - s " 



Xn{s)=X\ , n 
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so that Xn = 1 on /„. We use the pseudodifferential calculus and the notation 
defined in Section [3l taking the cut-off parameter r equal to r n (from condition 
(EIIR), as explained before the statement of theorem II .3|) . 

To simplify the notation, we will write ipj = ipji s j = S J i arL d a = a n - 
Thanks to Lemmas 13.101 and to the "Egorov property" Corollary 13.81 we havqf 

n n / * 

^2xn(sj) 2 \(^,a^ 3 )\ 2 = ^ I (t/>j , Op Gn (axn)i>j) 1 2 + nr 3 ( 

" / 1 

= I] |(V>,, Op G „(a T Xn)^)| 2 + nr 3 — — 

+o(t 2 ^ t \ J2 I \a{x^)\ 2 dv x {u:)xn{s) 2 d^s) 

+0{T 2 N 2 q r )\{x G D n ,p(x) <r + 2}\ f Xn(s) 2 dp(s) 
+0(T 2 ) J2 [ \a( s > N Hx,u)\ 2 dp x (u) Xn (s) 2 dn(s) 



where 



iefl n ' 



T-X _, iV-1 



= I^ aoff fc an d a ^N) = ^Y j q 2lsk aoa k . 



T 

a . 

T ^ N 

k=0 fc=0 



We use Lemma |3~31 to write 



•,T,, M|2 



3=0 

< ^ / |a T (a;,tj)| 2 d^(w)x„(s) 2 ^(s) 

ie6.D„ 

+?r X! / l« T (a ; ! w)| 2 ^(w) / Xn(s) 2 rff((s). 

xGD 7l ,p(x)<r 

We have seen in Section I2T21 that X^6,d„ / l aT ( :z: ; w)| 2 d^(w) = O ( ■£§ ). The same 

proof shows that J2 x eD n I \^ s ' N \x,u)\ 2 dv x (u) = O f -^ J . A major difference 
here with the usual Quantum Ergodicity (and with the special proof of 52) is that 
condition (EXP) is used already to show that the "remainder term" a^ s ' N ' of the 
Egorov theorem is small in the L 2 -norm. 

Remark 4.1. For s staying away from r/2, a slightly more careful proof would 
show that we only need to assume here that the spectrum of A is contained in 
{1} U [— 1, 1 — /3]. Hence our Remark 11.51 

Recall also, from the Kesten-McKay law (Section [SJ Corollary I5.2p that 

spfhsj /*.(•>■*<.) -0(1). 



°To prove the extended Theorem 1 1,71 we also need Lemma l3.11l 
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We obtain finally 

1 / 1 \ / /V 2 T 2 

E l^,«^>| 2 = ^ 3 ( T7 -r^ )+0' ■ 



N(I n ,G n ) s ^J^- n rj " \S n (rS n 

T 2 \ ( i 



and if we choose the sequences r = r n and S n as explained in section [5] 

Sj tin 

As the left-hand side of the equality does not depend on T and N, we take the limit 
N — > oo and then T — > oo to get 



lh ^ at(t r \ E |(^,a^-)| 2 =0. 



5. Kesten-McKay law for sequences of graphs satisfying (EIIR) 

In this section we give an alternative proof of the Kesten-McKay law |19L I26j , 
which gives the spectral density for large regular graphs satisfying (EIIR) and is 
analogous to the Weyl law for the spectral density of the laplacian on Riemannian 
manifolds. Note that we consider the density of eigenvalues in intervals that are 
allowed to shrink as n — > +00. 

In the definition of Op G , we take r = r n such that r n +2 < R n and q Tn a n — > 
(where R n and a n are the quantities occurring in (EIIR)jj. We also assume that 
there exists an integer M such that 

0. 



rM+1 M-3 
u n 1 n 

If \n is the function defined in (12. If) with sq e (0, r), this ensures that 

X„(s) 2 d//(s) > rl(5 n r n )- M . 


Theorem 5.1. Assume (EIIR). Let x — Xn be a smooth function satisfying 

fell < C k 5-\ 
with S n r n — > +oo, such that M+1 M _ 3 — > for some M . 



Then we have 



when n — > +00. 



1 

-^XniSj) 2 ~ / Xn{sfdll{s) 



I] 



^We can take for example r n = min < R n — 2, —(1 — e) ° gQ " I f for any < t < 1. 
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Proof. 

n 

Y,Xn( Sj ) 2 = Tr (Op G Jx«) 2 ) + 0(nr»(* n r„)-°°) 

3=1 

= E E I^Gj^y;Xn)| 2 + 0(nr3(«5„r„) — ) 

p(x)>r n yEX 

+ o{ q r ~) y, £!**(*, w;x»)i a 

p(x)<r n y£X 

= E Ei i ^<^«)i 2 + °( nr ™( <5 « r «^ 00 ) 

+ (0(«fn)-l) ^ £|i^(z,2/;X„)| 2 

p(x)<r n y£X 

= n f Xn(s) 2 d^( S ) + 0{nrl(5 n r n )-°°) 



+ (0(« r ») - 1) |{.t G D B ,p(x) < r„}| / X ™(s) 2 dM(s). 

Jo 



Thus, we get the desired result if 

X n {s)dn(s) > r£(d> n 7V, 



D 
Corollary 5.2. Under the assumptions of Theorem \l.S\ we have 

N(I n , G n )~n [° " d M («) - 27^ T1 |c(s )|- 2 

JSO— in 

where \c{s$)\ is the density of the Plancherel measure at s . 

6. Quantitative statement 

In this section we will give explicit upper bounds on the rate of convergence, 
first in terms of the parameters R n and a n associated with the sequence of graphs 
(G n ) in condition (EIIR), then depending only on n for sequences of random graphs. 
These results are certainly not optimal because some of our inequalities were written 
in a non optimal way. 

In the general case so € (0, r), we have 

Lemma 6.1. If S n = r~ 1+e for some < e < 1, then we have 

j^ e |«r>,«„<>>f=o (■•-'» +r r '-"«•■)• 

s ( " } ei n 
where we can take r n = min {i?„ — 2, —(1 — e') log„(a n )} , for any < e' < 1. 



QUANTUM ERGODICITY ON LARGE REGULAR GRAPHS 25 

Proof. According to the proof of section 0J we have 



wb^S ^'•^' ) ^ rl0 {l^F) +0 



N 2 T 2 






Take A = rj and T = rj such that i = 2L = N Jf = r,7 . For every 

M > 0, we have O ( Sn(r r } Sn) «) = O (r3- M <5 r 7 (1+M) ) = O (r^ (1+M)e ) and this 
term can be made negligible in comparison with the other terms by taking M 
sufficiently large. Finally N 2 T 2 q 7 ' n a n — r n q r "a n . D 



Remark 6.2. Here we kept the spectral gap /3 fixed, but we see that this could 

2/9 

be relaxed to (3 >• r n 

For sequences of random graphs, we have 

Lemma 6.3. Let 8 > 1/2, e > 0, and S n = (log q (n 1 ~ s )) 1 ~ t . If G„ is chosen 
uniformly at random among the (q + l)-regular graphs with n vertices, we have 

with overwhelming probability. 

Proof. Take i? n and a n as in example [T] Let 1/2 < 8 < 1, then i? n = 
(1 — 8) log (n) and a n = 80(1 — <5) log„(n)ri _l5 . In this case, we take 

r„ = min(i?„ - 2, -(1 - e ')i2^H 1 = i?„ _ 2 , 

and apply lemma [6~T1 D 

7. Proof of Theorem 11.71 

Most steps of the proof carry over to arbitrary a £ S® (X) . Actually, all that 
needs modifying is the treatment of the expression 

2 



(x, lu, s)dv x (u>)dfx(a) 



xeD n J k=o 

that is used in M2.2I (for s = 0) and Section 0] (for s close to s ). Equation (|7.1[) is 
also 

1 „Af-l fc 

(7.2) ^jE l^2^2Q 2lsJ a ^ : '(x,u,s)a(x,u},s)dv x (u>)d^(s) 

■ " - '- =0 3=0 

1 „JV-1 fc 



A 2 

xGD„ fc=0 i=l 



In 92~21 the integral 

> / a o er^ (x,a;)a(a;, uS)dv x (ui) 



xe_D„ 
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was rewritten as J2xeD Sjd(x)a(x) using the fact that a did not depend on uj - 
thus establishing a link between the shift a and the laplacian. We need to adapt 
that argument to the case when a{x,uj) depends on the first D coordinates of the 
half geodesic [cc,w) = (x,X\,x%,. . .). 

7.1. Proof when D = 2. When D — 2, a(x,uj) = a(x,Xi), so that a is a 



B 



if e is an element of B, we shall denote by o(e) € V n its origin, t(e) £ V n its 
terminus, and e *E B the reversed bond. 
One sees that 

} / ao a-'fx, uj)a(x,u))dv x (uj) = } M^a(e)a(e) 

J q + 1 

xeD n J H eeB 

where M* is a bistochastic matrix indexed by B, defined by 

Af J (e,e') = - 

if o(e') = t(e) and e' ^ e; and M"(e, e') = otherwise. This is (up to normalization) 
the matrix appearing in §3 of [30j . It is the (normalized) adjacency matrix of the 
g-regular directed graph, whose vertices are the directed bonds of G, and where 
we draw an edge between two bonds if they are consecutive without allowing back- 
tracking. What we need is an explicit relation between the spectrum of M* and 
the spectrum of the discrete laplacian on G, in other words, of the matrix A. The 
relation between the eigenvalues is formula (44) in |30j . but since we also need 
relations between the eigenf unctions, we shall be more explicit below. We did not 
write all the detailed calculations because they are lengthy but basic. We assume 
these relations must already be known but did not find any reference. 

(o) Both A and M" have 1 in their spectrum, corresponding to the constant 
eigenfunction. The matrix A has —1 in its spectrum iff the graph G is 
bi-partite, in which case M" also trivially has —1 in its spectrum. 

(i) each eigenvalue A 7^ ±1 of A gives rise to the two eigenvalues 

2 



of M"; 
(ii) in addition, M" admits the eigenvalue 1/q with multiplicity b :— \E n \ — 
\V n \ + 1 (the rank of the fundamental group of G); and the eigenvalue 
— 1/q with multiplicity 6—1 if —1 is not an eigenvalue of A, or b if — 1 is 
an eigenvalue of ArH 
In particular, the eigenvalue 1 of M" has multiplicity 1. The tempered spectrum 
of A corresponds to eigenvalues of Af" of modulus 1/y/q; the untempered spectrum 
of A contained in [—1,1-/3] gives rise to real eigenvalues of M" contained in 
[-1,1-/3'] with 

i-r 



(?+i)(i-.s- v /(i-« 2 - 5 s : 



°Note that B has cardinality n(q + 1) if G has n vertices and is (q + l)-regular. 
Or, equivalcntly, if the graph is bi-partitc. 
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Since M* is not normal, the knowledge of its spectrum is not sufBcient to control 
the growth of M* k in a precise manner (we need a bound that is independent of the 
size of the matrix, in other words, independent of n). Below, we describe explicitly 
the eigenvectors of M" in terms of those of A; these eigenvectors do not form an 
orthogonal family but this is compensated by the fact that one can compute their 
scalar products explicitly. 

(i) an eigenfunction <p of A for the eigenvalue A ^ ±1 gives rise to the two 
eigenfunctions of M' , 

/i(e) = 0(t(e)) - ei0(o(e)); / 2 (e) = #t(e)) - e 2 <f>(o(e)), 

where ei, £2 are the two roots of qe 2 — (q + l)Ae + 1 = (in what follows 
we index them so that |ej.| < | e2 1 ) - Special attention has to be paid to the 
case A = ±-¥y, for which ei = e 2 (see below), 
(ii) the eigenvalues ±l/g of M" correspond, respectively, to odd and even 12 ! 
solutions of J2ote)=x /( e ) = (f° r every vertex x). For the eigenvalue 1/q, 
an explicit basis of eigenfunctions is indexed by generators of the funda- 
mental group, (71,..., 7;,) : every closed circuit 7 made of consecutive 
edges (ei, . . . , e^) gives rise to an odd eigenfunction 

k 

f-y = J2 ^ - s *i ■ 

If G is bi-partite, then all circuits have even length and we have an explicit 
basis of even eigenfunctions for the eigenvalue —1/q, again indexed by 
generators of the fundamental group : 

k 

ff7=E(-l)*(^+^) 

i=i 

if 7 is a closed circuit made of consecutive edges (ei, . . . , e^). If G is not 
bi-partite, there are closed circuits of odd lengths, in which case g 7 is 
not an eigenfunction of M*. Nevertheless, if 7,7' are two circuits of odd 
lengths, g 1 — g 7 > is now an eigenfunction of M* for the eigenvalue —1/q. 

The eigenfunctions of the family (ii) are automatically orthogonal to those of the 
family (i). In (i), eigenfunctions of M* stemming from different eigenvalues A of A 
are orthogonal; however, the two eigenfunctions fi, f 2 stemming from the same A 
are not orthogonal. 

To evaluate the norm of a matrix, it is safer to work in an orthogonal basis, 
and thus we shall consider, instead of a pair (/i,/2), the pair 

h(e) = ^(t(e)) - e x 4>(o(e)), f 2 (e) = 0(t(e)) - /^(o(e)) 

which can be checked to be orthogonal for 

eiA- 1 



/' 



ei-A 



12 



Odd means /(e) = —/(e) and even means /(e) = /(e), for every bond e. 
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Iii the plane generated by I ^tt , nth ) , M$ has matrix 

1 to J \\\fi\\ ' u/aiiy 

1/gei * 
l/qe 2 

where * is a number that can be calculated explicitly in terms of ei, £2 and A, and 
which is uniformly bounded (since the norm of M", anyway, is bounded indepen- 
dently of n). 

This discussion is also valid for A = ±—7-?, a special case where t\ = €2 = ±-4=. 

To summarize, the spectrum of M* is contained in [— 1, 1 — (3'} U {1} (resp. 
[-1 -f/J', 1 - 0'] U {1}) if the spectrum of A is contained in [-1, 1 - 0\ U {1} (resp. 
[— 1 + /3, 1 — /3] U {1}). We can find an orthonormal basis of L 2 (C B ) in which M" is 
block diagonal, each diagonal block being an upper triangular matrix of size < 2, 
and the non-diagonal coefficients are uniformly bounded. 

This implies that the operator 

N-l k 

Jp EE^ 1 ^' = ■^(q 2ls M i -I)- 2 (g 2 -( JV + 1 )M»( iV+1 ) - q 2ls M* - N(q 2i8 M* - I) 

k=Q j=0 

has norm O (-^) on the orthogonal of the constant function (for any real s if the 
spectrum of M" is contained in [—1 + /3', 1 — /?'] U {1}, or for q 2ls away from —1 if 
the spectrum of M' is contained in [— 1, 1 — f3'] U {1}). This tells us that (|T.2[) . and 
hence ([73]). is O(i). 

7.2. Reduction to the case D = 2. Let us now consider Theorem II .71 in the 
case of an operator whose kernel on the tree satisfies K (x, y) 7^ => dx(x, y) = D. 
Theorem 11.71 proven in the case D — 2, can be applied to the (q + l)g- D_1 -regular 
graph with vertex set V n and adjacency matrix (q + l)q D ~ 1 SD, with the notation 
of (|2.2p . This implies Theorem 1 1.71 in the general case. 
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